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An approximate equation was derived to predict the mass of insulation required to 
limit the maximum temperature reached by an insulated structure subjected to a transient 
heating pulse. In the course of the derivation two figures of merit were identified. One figure 
of merit correlates to the effectiveness of the heat capacity of the underlying structural 
material in reducing the amount of required insulation. The second figure of merit provides 
an indicator of the mass efficiency of the insulator material. An iterative, one dimensional 
finite element analysis was used to size the external insulation required to protect the 
structure at a single location on the Space Shuttle Orbiter and a reusable launch vehicle. 
Required insulation masses were calculated for a range of different materials for both 
structure and insulator. The required insulation masses calculated using the approximate 
equation were shown to typically agree with finite element results within 10 to 20 percent 
over the range of parameters studied. Finite element results closely followed the trends 
indicated by both figures of merit. 


Nomenclature 

c pe Effective insulator specific heat capacity 

c ps Structural specific heat capacity 

d e Insulator thickness 

d s Structural thickness 

k e Effective insulator thermal conductivity 

m Mass per unit area 

m e Insulation mass per unit area 

m s Structural mass per unit area 

m op t Minimum total mass 

m so Structural mass at which total mass is minimum 

Pavg Average ambient pressure during equivalent square surface temperature pulse 
T Temperature 

t Time 

Th Applied surface temperature rise of insulator 
th Duration of heating pulse 

Ti Initial temperature 

T m Maximum structural temperature rise 

T ce Temperature to use for calculating effective insulation specific heat capacity 

T cs Temperature to use for calculating effective structural specific heat capacity 

Tke Temperature to use for calculating effective insulation conductivity 

Symbols 

/3 S Figure of merit for thermal performance of structure 
7 Ratio of insulator to structural heat capacitance 
K e Figure of merit for thermal performance of insulator 
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p e Effective insulator density 

p s Structural density 

Th Non-dimensionalized duration of heat pulse 

I. Introduction 

A EROSPACE vehicles in hypersonic flight generally require some form of thermal protection system to 
accommodate the aerodynamic heating on the vehicle surface. One approach is to cover the exterior 
surface with a thermal insulator, such as the reusable ceramic tiles and blankets used on the Space Shuttle 
Orbiter. Typically, the vehicle structure would be designed for structural considerations. The insulation 
would then be sized to maintain the structure below the maximum allowable temperature for all design 
environments and loads. However, the amount of required insulation is a function of the mass and properties 
of the underlying structure. Also, it is not clear what combination of insulation properties would lead to the 
lightest insulation for the transient thermal surface heating. 

In this paper an approximate solution ’ ’ for maximum structural temperature of an insulated structure 
is used to explore parameters that drive the combined mass of structure and insulator. Equations are 
developed for the minimum insulator mass required to maintain the inner, unheated surface below a specified 
temperature. A figure of merit which correlates to the effectiveness of the heat capacity of the underlying 
structural material in reducing the amount of required insulation was developed. A second figure of merit 
was identified for the combination of insulator thermal properties that minimizes the mass of an insulator 
required to limit the maximum temperature of an underlying structure subjected to a transient heating 
pulse. Insulation was sized for a number of surface heating histories, insulators, and structural materials by 
iterating a one-dimensional, finite element analysis and by using the simple equations derived in this paper. 
Finite element results were correlated to the two figures of merit developed in this paper. 


II. Problem Definition 


The simplified thermal problem that 
was analyzed by Blosser ,2 is illustrated Heated Surface 

in Figure 1. A thermal insulator with 0 t 
thickness, d e , density, p e , specific heat 
capacity, c pe , and thermal conductivity, 
k e covers a structure (inner face sheet) > ^ 
of thickness, d s , density, p s , and specific 
heat capacity, c ps . The inner surface of 
the structure is assumed to be perfectly 
insulated to simplify the mathematics of 
the problem and because it is a com- 
monly used conservative assumption for 
sizing thermal protection systems. To ci ] 
further simplify the problem, the struc- 
ture is treated as a lumped heat capac- 
itance. For this solution the material 
properties are assumed to be constant, so 

effective, averaged properties would have to be used to approximate more complex material behavior. 

A simple transient heat pulse was assumed. Initially the insulator and structure are assumed to be at 
a uniform temperature of 0. The outer surface of the insulator is assumed to instantaneously rise to a 
temperature Th at t = 0 and maintain that temperature until t = th at which time it instantaneously returns 
to 0. 
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Figure 1. Illustration of simplified problem 


Time 


III. Maximum Structural Temperature 

An approximate solution developed for the maximum structural temperature of the simplified problem 
illustrated in Figure 1 is given by: 
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where 


( 2 ) 


and 
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Th = 


Pe^pe^e 

PsCps^s 

k e 


PeCpe 


dl 


th 


( 3 ) 


IV. Mass of Insulated Structure 

For aerospace vehicles reducing mass is of utmost importance. Therefore, it is helpful to develop further 
insight into the interplay between the mass of structure and the mass of insulator required to limit the 
structural temperature. The mass per unit area of this simplified, insulated structure can be written simply 
as: 

m = m e + m s = p e d e + p s d s (4) 


The relationship between d e and d s that is required to limit the maximum structural temperature to 
a specified value is given by Equation (1). By substituting Equations (2) and (3) into Equation (1) and 
rearranging terms, the following expression can be obtained: 


d s di = 


( k e t h ) 2 


2 (Ps^ps) (PeCpe) ^ lu{\ 


( 5 ) 


The equation can be further manipulated to give an explicit expression for d e . This provides a very 
useful equation to calculate the insulator thickness required to limit the maximum structural temperature 
to a specified value. 


d e 


( k e t h ) 2 


y 2(/9 s Cp S d s ) (p e Cp e ) ^ ^77.(1 rp^ 


( 6 ) 


Substituting Equation (6) into Equation (4) produces an expression for mass per unit area as a function 
of structural thickness: 



( 7 ) 


The first term in Equation (7) represents the mass of insulator required to limit the structure to the 
specified maximum temperature and the second term represents structural mass. Equation (7) can be 
written as a function of m s . 


m = 



( 8 ) 


IV. A. Material Properties and Insulator Mass 

The first term in Equation (8) represents the mass of insulator required to limit the structure to the 
specified maximum temperature. Inspection of this equation illustrates the effect of key parameters on 
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required insulator mass. 



( 9 ) 


Table 1. Simplified surface temperature histo- 


From Equation (9) it is obvious that a longer heat pulse (increased th) will result in a larger required 
insulator mass. Similarly, a larger structural mass will provide more structural heat sink capacity and require 
less insulator mass. For a given structural mass and a particular transient heating pulse, required insulation 
mass will be determined by the thermal properties of the insulator and structure. 

The specific heat capacities of insulators and structural ma- 
terials vary significantly with temperature and the thermal 
conductivity of low density insulators is a function of both 
temperature and pressure. Therefore, it is not immediately 
obvious which material property values to use in the preced- 
ing equations to obtain results that are meaningful and useful. 

A method was previously devised to calculate an equivalent 
square heating pulse that approximates a surface temperature 
history on a vehicle during atmospheric entry, resulting in val- 
ues of th , Th and P avg ■ These values are shown in Table 1 for 
the seven heating histories previously studied. ' The first four 
temperature histories were measured at four different locations 
on the windward surface of a Shuttle Orbiter during a single 

atmospheric entry flight. *’ It is interesting to note the variation in the equivalent square heating pulses that 
occur over just the windward surface in a single flight. The last three heating histories were predicted for a 
single point on the windward surface of three different proposed reusable launch vehicles. ,l1 Also shown is 
the temperature at which the effective insulator thermal conductivity should be calculated. 


Profile 

th 

S 

T h 

°C 

P aV9 

Pa 

Tke 

K 

BP9470 

1273.6 

838.6 

98.0 

792.0 

BP9678 

1313.8 

799.9 

133.0 

768.8 

BP9591 

1206.6 

804.0 

99.7 

771.3 

BP9489 

1191.2 

576.0 

59.1 

634.5 

ATSpA 

1914.4 

871.2 

1756.6 

811.6 

RLVpA 

2042.1 

861.8 

1025.5 

806.0 

RLV3c 

1919.6 

727.0 

1364.6 

725.1 


The = Ti + 0.6 Th, where T, = 60°F, 289A' 


(10) 


Each of the relevant thermal properties was examined independently to devise a method to calculate an 
equivalent constant value corresponding to a particular transient surface temperature history. The equivalent 
constant values for a particular surface heating history can be calculated as follows: 

c ps = c ps (T cs ), where T cs =Ti + (11) 

Cpe — (T c e J , where T ce — Tj -f- T m: and ( \ ‘2j 

k e = k e {T kei P avg ). (13) 


Inspection of Equations (8) and (9) reveals that the thermal properties of the structure and insulator 
can each be grouped separately. These grouped properties have the potential to be useful for comparing the 
thermal effectiveness of both the structural and insulator materials. 

The specific heat capacity and maximum allowable temperature of the structure directly affect the re- 
quired insulation mass. However, the thermal properties of the structure are coupled to the amplitude of the 
equivalent square heating pulse, T k - The following expression is a candidate figure of merit for the thermal 
effectiveness of the structural material: 



A larger value of this parameter will lead to a smaller required mass of insulation. Therefore a structure 
with a high specific heat capacity and high maximum temperature limit will tend to require less insulation. 
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For steady state heat conduction, minimizing the product of insulator density and thermal conductivity, 
p e k e , minimizes the mass of required insulator. However, a similar figure of merit for minimum mass of 
insulation subjected to transient heating has not been available. The grouping of insulator properties in 
Equation (9) suggests a candidate figure of merit for minimum mass insulation for a transient heating pulse 
(or at least for the simplified heat pulse in this derivation): 


K e = 



(15) 


A smaller value of this parameter will lead to a smaller required mass of insulation. Therefore, a low 
mass insulator for the transient heat pulse should have a combination of low p e k e and high c pe . 


IV. B. Minimum Mass of Insulated Structure 


The mass of the insulated structure, calculated using Equation (8), is obviously a function of the mass 
of the structural material. Increasing the amount of structural mass linearly adds to the total mass, but 
results in a corresponding decrease in the mass of required insulation. Therefore, in the absence of any 
overriding structural requirements, there should be a structural mass for which the total mass is minimum. 
The structural mass for which the total mass is minimum can be readily calculated by taking the first 
derivative of Equation (8) with respect to m s and setting it equal to zero. The result can be solved for m s 
to produce: 


= 54 ( ^y/th 


'-'PS 


te ) 2 

{~ lrl ( 1 “ Tf)) 


(16) 


Equation (16) is potentially useful for sizing insulated structure. The structure could first be sized to 
carry the design structural loads. The resulting equivalent structural mass can then be compared to that 
calculated using Equation (16). If the result of Equation (16) is greater, then increasing the structural mass 
to that value will result in the lowest combined mass of structure and insulator. However, if the result of 
Equation (16) is less, then adding additional structural mass will only increase the combined mass. 

Substituting Equation (16) back into Equation (8) produces the following expression for the minimum 
mass of an insulated structure: 


rriopt = (3 + 1) m so 


(17) 


TMopt 



(18) 


The first term in parentheses in Equation 17 indicates the relative contributions of the insulation and 
structure, respectively, to the minimum total mass. For minimum total mass, Equation 17 predicts that the 
mass of the insulation will be three times the mass of the structure. Equation (18) defines the lower bound 
for the mass of an insulated structure subjected to a transient heat pulse. The magnitude of the minimum 
total mass is solely a function of the duration of the equivalent square heating pulse, th , and the figures of 
merit defined in Equations (14) and (15). 


V. Analytical Predictions and Numerical Simulations 

Several numerical studies were performed to assess how the simple approximate solutions derived in 
this paper compare to more accurate numerical simulations of earth atmospheric entry heating on reusable 
launch vehicles. The same finite element model and heating histories'’ were used for all of these numerical 
studies. The objectives of the numerical studies were: 1) to assess the accuracy of using Equations 6 and 
9 to size insulation for a realistic transient heating pulse, 2) to determine the applicability of the figures of 
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merit defined by Equations 14 and 15 to results of atmospheric entry simulations, and 3) to investigate the 
interplay between the mass of the structure and the corresponding mass of insulation required to limit the 
maximum structural temperature. 

V.A. Finite Element Model 

A one-dimensional finite element model was de- 
veloped using the DOLFIN finite element library 
for the Python programming language. Using the 
DOLFIN library, linear one dimensional elements 
were used to discretize the spatial dimension and 
an implicit Crank-Nickelson time marching scheme 
was used to solve the weak formulation of the dif- 
fusion equation. The model consisted of a layer of 
insulator material in perfect contact with a layer 
of structural material similar to the configuration 
shown in Figure 1, except that the structure in the 
finite element model was not treated as a lumped 
heat capacitance. For the results shown in this pa- 
per, the model consisted of 50 elements through the 
thickness of the insulator and 2 elements through 
the thickness of the structure. The boundary condi- 
tion on the outer surface of the insulator consisted 
of an imposed surface temperature that could be varied arbitrarily with time and updated at each time step 
of the solution. The inner surface of the structure was adiabatic. The material properties could be arbitrary 
functions of temperature and ambient pressure. Material property values for each finite element were up- 
dated at each time step of the analysis. Property values could be different for each element, but did not vary 
spatially within an element. Time steps between 1 and 5 seconds were used to calculate the results presented 
in this paper. The finite element simulation was iterated to determine the insulation thickness required to 
limit the structural temperature rise to the specified value. Each transient simulation was continued until 
the structural temperature had reached its maximum value. 


Table 2. Thermal Properties of Four Structural Materials 


Ti + T m 

c ps (T cs ) 

/3„( BP 7940) 

/3 s (ATSpA) 

K(°F ) 

~WgK ( IbmR ) 

J ( Btu \ 
kg K ^ IbmR ) 

J ( Btu \ 
kgK IbmR J 

Aluminum 2024 p s 

=2803 ^ (0.10 

422(300) 

450(350) 

900(0.215) 

910(0.217) 

27(0.006) 

41(0.010) 

25(0.006) 

38(0.009) 

Graphite/epoxy p s 

= 1576 (0.057 l -i^) 

366(200) 

422(300) 

913(0.218) 

984(0.235) 

9(0.002) 

30(0.007) 

8(0.002) 

27(0.006) 

Beryllium Aluminum 

p 3 =2098 

(0.076 ^J) 

505(450) 

561(550) 

1690(0.404) 

1735(0.415) 

151(0.036) 

267(0.064) 

138(0.033) 

244(0.058) 

Titanium 6A1-4V p 

s =4437 ^ (0.16 ^J) 

533(500) 

700(800) 

574(0.137) 

594(0.142) 

68(0.016) 

270(0.064) 

62(0.015) 

242(0.058) 


V.B. Thermal Properties of Structural Material 


I Ibm 

ft 1 


Obviously, raising the structural tempera- 
ture limit will decrease the amount of required 
insulation. Less commonly considered is the 
effect of the structural specific heat capacity 
on the required insulation. Equation 14 de- 
fines a candidate figure of merit that combines 
the effects of structural temperature limit and 
specific heat capacity on the amount of re- 
quired insulation. 

A numerical study was performed to as- 
sess how well the results of a series of full 
numerical simulations would correlate to the 
parameter /3 S . To compare the effects of the 
structural thermal properties, it was neces- 
sary to fix all other parameters in the problem 
and vary only the choice of structural mate- 
rial. LI-900 tile was chosen as the insulator 
and the structural mass per unit area was 
fixed at 6.10 ^ (1.25 jjt )- Two different 
heating histories were considered, but results 

can only be directly compared for a particular heating history. Table 2 shows thermal properties for four 
structural materials chosen for the study. The temperature dependent properties for these materials are 



Figure 2. Effect of structural thermal properties on required insu- 
lator mass 


6 of 10 


American Institute of Aeronautics and Astronautics 


tabulated in another paper. The maximum allowable temperature for each material is application depen- 
dent, so an approximate upper and lower bound was chosen for each material. For each limiting temperature 
the structural specific heat capacity (Equation 11) and j3 s for each the two heating histories are shown in 
Table 2. 

Iterative finite element analysis was used to size insulation for 16 cases (4 structural materials, 2 temper- 
ature limits, and 2 heating histories). The resulting insulation masses are shown in Figure 2 and indicated 
by the square symbols. For the ATSpA heating history, each particular structural material and temperature 
limit is identified in the figure. The square symbols are not labeled for the BP7490 heating history to avoid 
cluttering the figure. The approximate insulation masses, calculated using Equation 9, are indicated by 
circles corresponding to each finite element solution. For the Space Shuttle heating history, BP7490, half of 
the approximate solutions are within 10 percent of the finite element solutions and all of the approximate 
solutions are within 25 percent. For the reusable launch vehicle heating history, ATSpA, half of the approx- 
imate solutions are within 10 percent of the finite element solutions and all of the approximate solutions 
are within 20 percent. Because the value of c pe used in Equation 9 is a function of the maximum structural 
temperature limit (Equation 12), the approximate insulation masses are not solely a function of /3 S . The 
solid lines in Figure 2 represent the approximate solution using a value of c pe that is an average of the 8 
cases for each temperature history. The close agreement between the solid lines and circles indicates that, 
for LI-900 insulation, the variation of c pe with T m has little effect on the mass of required insulation. The 
solid lines clearly indicate the effect of /3 S on the required insulation mass predicted by the approximate so- 
lution. Remarkably, the finite element solutions closely follow this trend, even with four different structural 
materials. Therefore, the parameter f3 s appears to be a useful indicator of the thermal effectiveness of a 
structural material. 

V.C. Thermal Properties of Insulator 

To better understand the significance of the pa- 
rameter defined in Equation 15 it is helpful to con- 
sider a range of candidate insulator materials that 
could be used for thermal protection on a reusable 
launch vehicle. A numerical study was performed 
to investigate the relative mass efficiency of a range 
of insulators for two different surface temperature 
histories (BP7490 and ATSpA) from Table 1. Five 
rigid ceramic tile materials (LI-900, LI-2200, FRCI- 
12, AETB-8, AETB-16) and four flexible insulations 
(Safhl-3, Saffil-6, QHber-3, and Qfiber-6) were con- 
sidered. The temperature and pressure dependent 
thermal properties of the three Space Shuttle Or- 
biter tile materials ((LI-900, LI-2200, and FRCI-12) 
were obtained from published tables. The thermal 
properties for AETB-8 and AETB-16 were calcu- 
lated using the model developed by Daryabeigf 1 and 
the flexible insulation properties were calculated us- 
ing the models also developed by Daryabeigi. For 
this study, aluminum with a maximum temperature 
limit of 450K (350°F) was chosen as the underlying 
structure. The structural mass per unit area was 
fixed at 6.10 % (1.25 l jf). 

Table 3 shows the effective c pe (Equation 12) and 
k e (Equation 13) of each insulation for the two sur- 
face heating histories. The insulation figure of merit, 
calculated using Equation 15, is shown for each set 
of material properties. As shown in Table 3, the insulation effective material properties and figure of merit, 
K e , for a given insulation can be significantly different for different surface temperature histories. Therefore, 
insulations should be directly compared only for a given surface temperature history, structural material, 
and structural mass. 


Table 3. Insulation figure of merit values for aerospace in- 
sulations 


Temp. History 

Cpe 

fee K e = 4^ 

y/ c pe 


W 

J 

kg 2 


mK 

kg K 


LI-900 / 

=144 J (9 J) 

BP7490 

949 

0.03399 

0.1589 

ATSpA 

949 

0.06060 

0.2832 

LI-2200 / 

=S52 J (22 % 

f) 

BP7490 

949 

0.06046 

0.6907 

ATSpA 

949 

0.08326 

0.9512 

FRCI-12 

p s =192 ^ (12 ^ 

T ) 

BP7490 

949 

0.04039 

0.2517 

ATSpA 

949 

0.07212 

0.4494 

AETB-8 

p„ =128 J ( 7 8 773 

) 

BP7490 

975 

0.05702 

0.2338 

ATSpA 

975 

0.08728 

0.3578 

AETB-16 

Ps =256 J (16 J) 

BP7490 

975 

0.05927 

0.4860 

ATSpA 

975 

0.08975 

0.7360 

Saffil-3 

II 

00 

5 < 3 


BP7490 

975 

0.04352 

0.0669 

ATSpA 

975 

0.07302 

0.1123 

Saffil-6 

p s =96 

1 J (6 J) 


BP7490 

975 

0.02420 

0.0744 

ATSpA 

975 

0.05489 

0.1688 


Q fiber-3 p s =48 J (3 (J) 


BP7490 

ATSpA 

932 

932 

0.04106 

0.07328 

0.0646 

0.1152 

Qfiber-6 

p s =96 

5( 6 t? 

) 

BP7490 

932 

0.02158 

0.0679 

ATSpA 

932 

0.04975 

0.1564 
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The iterative finite element analysis, previously described, was used to calculate the required insulation 
mass for each insulation and surface heating history listed in Table 3. Figure 3 shows the required insulation 
masses versus the insulation figure of merit, n e . Each symbol shape represents the finite element solution 
for each of the insulations represented in Table 3 and the solid line represents the approximate solution 
calculated using Equation 9. The red line and symbols correspond to the BP7490 temperature history and 
the black line and symbols correspond to the ATSpA temperature history. For each temperature history, the 
finite element results for this wide range of insulators clearly follow the trend predicted by the approximate 
solution. 

For both surface temperature histories the 
flexible insulation materials are more mass ef- 
ficient than the rigid tile materials. For the 
Shuttle Orbiter surface temperature history 
(BP7490), the bulk of the heating takes place 
at such a low pressure that gas conduction 
is not significant and radiation dominates the 
heat transfer through the insulation. There- 
fore the density of the flexible insulation (af- 
fects the length of gas conduction path) has 
little effect on the thermal performance, so 
that all four flexible insulations have nearly 
the same mass efficiency. Conversely, the bulk 
of the heating for the reusable launch vehicle 
surface temperature history (ATSpA) occurs 
at a higher pressure at which gas conduction 
through the insulation is significant. There- 
fore, the lower density flexible insulations are 
more mass efficient than the higher density 
flexible insulations. For both surface temper- 
ature histories, lower density ceramic tiles tend to be more mass efficient than higher density tiles because 
heat transfer due to solid conduction becomes more significant with density. However, the denser LI-900 tile 
is more mass efficient than the AETB-8 tile. The required insulation masses predicted by the approximate 
equation are within 15 percent of the finite element solutions. 

V.D. Interplay Between Mass of Insulator and Structure 

An additional numerical study was per- 
formed to gain a better understanding of the 
interplay between the mass of the structure 
and the mass of insulator required to limit 
the structural temperature. LI-900 ceramic 
tile was used as the insulator. Aluminum, 
with a maximum temperature limit of 450K 
(350° -F), was chosen as the underlying struc- 
ture. The structural mass per unit area was 
varied between 0.49 and 14.6 (0.1 and 3.0 

y|J ) . For each of the two different surface tem- 
perature histories (BP7490 and ATSpA) from 
Table 1, iterative finite element analyses were 
performed for thirteen different values of struc- 
tural mass per unit area to calculate the re- 
quired insulation mass. 

Results of this numerical study are shown 
in Figure 4. The red symbols and lines corre- 
spond to the BP7490 heating history and the 
black to the ATSpA heating history. The blue 
line is the mass of the structure, which was the same for both heating histories. Required insulation masses 
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Figure 4. Mass of insulated structure as a function of structural 
mass 
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Figure 3. Required insulation mass for several insulators 
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calculated using the finite element analyses are represented by the solid circles. Dashed lines represent the 
corresponding approximate insulation masses calculated using Equation 9 with properties defined by Equa- 
tions 11-13. The solid squares represent the sum of the finite element solutions for insulation mass and the 
structural mass. The solid lines represent the sum of the approximate solution for insulation mass and the 
structural mass. The solid triangles represent the minimum total mass calculated using Equations 16 and 
18. 

Comparison of the required insulation masses calculated using the finite element analysis and using 
Equation 9 (solid circles vs dashed lines in Figure 4) shows remarkably good agreement. Over the middle 
range of structural masses studied the approximate solution is within 15 percent of the finite element solution 
for both heating histories. However, for very low values of structural mass the approximate solution begins 
to sharply diverge from the finite element solution. Also, at the high end of the range of structural masses, 
the approximate solution for the BP7490 heating history is beginning to diverge from the finite element 
solution. This indicates that the approximate solution may become much less accurate when the ratio of 
insulation mass to structural mass becomes very high or very low. The total mass is obtained by simply 
adding the structural mass to the insulation mass, so the accuracy comparison between the finite element 
and approximate solutions is similar to that for the required insulation masses. 

As the structural mass increases the required insulation mass decreases, so there is the potential for 
a combination of structure and insulator that produces the minimum total mass. For the finite element 
results the total mass minimum occurs at a structural mass value of about 1 (0.2 jft) for both heating 

histories. Equation 16 predicts a minimum at 2.1 ^ for the BP7490 heating history and 3.6 for the 
ATSpA heating history. Although Equation 16 does not accurately predict the structural mass at which the 
total mass is a minimum, the values of the minimum total mass predicted using Equation 18 are within 10 
percent of the finite element calculations for both heating histories. The three to one ratio of insulator mass 
to structural mass predicted to occur (Equation 17) at the minimum mass appears to be an artifact of error 
in the approximation. 

The total mass curves in Figure 4 are relatively flat over a significant range of structural masses. There- 
fore, the mass (and thickness) of the insulator can be reduced by increasing the structural mass, for a 
relatively modest total mass penalty. Thinner vehicle walls, insulator plus structure, could improve vehicle 
packaging efficiency, perhaps resulting in a net mass savings. 

VI. Summary 

Approximate analytical equations were developed for sizing the insulator thickness and mass required 
to maintain the insulated structural skin of hypersonic aerospace vehicles below a specified temperature. 
Manipulation of the equations revealed that the thermal properties of the insulator and of the structure 
could each be collected into a single term. These terms could be used as figures of merit to indicate the 
effect of the choice of insulator or structural material on the required insulator mass. An equation was also 
developed for the minimum total mass of an insulated structure sized to stay below the maximum structural 
temperature limit, ignoring other design considerations. 

A one-dimensional, finite element, transient thermal analysis was used to perform numerical studies for 
comparison with the results of approximate equations developed in this paper. The finite element analysis 
included the effects of temperature and pressure dependent material properties. Time-accurate surface 
temperature histories and ambient pressure histories were also incorporated into the numerical models. 

A numerical study was performed to compare the mass of insulation required to protect 4 different struc- 
tural materials. Two different maximum structural temperatures were considered for each structural mate- 
rial. The insulation material and the structural mass were kept the same, but calculations were performed for 
two different surface temperature histories. The approximate solution predicted required insulation masses 
within 10 percent of the finite element results for more than half of the cases and within 25 percent for all 
of the cases considered. The finite element results exhibited a close correlation to the figure of merit, f3 s , 
derived in this paper. 

A second numerical study was performed to compare the mass efficiency of 9 different aerospace insula- 
tors, five ceramic tile materials and four flexible insulations. The structural material, mass and maximum 
temperature limit were held constant while the required insulation masses were calculated for the same two 
surface temperature histories. For the insulations considered, the required insulation masses predicted by 
the approximate equation were within 15 percent of the finite element solutions. The finite element results 
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were closely correlated to the insulation efficiency figure of merit, n e , derived in this paper. 

A third numerical study was performed to investigate the effect of structural mass on the mass of re- 
quired insulation and to determine the minimum mass for an insulated structure designed with only thermal 
constraints. The structural material and maximum temperature, as well as the insulation material, were 
held constant and the required insulation masses were calculated for a range of structural masses and for the 
same two surface temperature histories. For very small or very large values of structural mass the approxi- 
mate solution diverges from the finite element solutions, however, away from these extremes the approximate 
solution was within 15 percent of the finite element solution. The approximate equation for minimum total 
mass was within 10 percent of the finite element solution, however, the approximate equation predicted the 
minimum total mass at a much larger value of the structural mass than the finite element solution. 

The approximate equations developed in this paper were shown to predict the results of much more 
complex finite element calculations with surprising accuracy for the wide range of cases considered in this 
paper. Achieving these accurate results requires carefully following the techniques developed for mapping 
the transient surface temperature histories to equivalent square temperature pulses and for calculating the 
effective property values to substitute into the approximate equations. These simple approximate equations 
have the potential to be very useful for the preliminary investigation of a wide range of design space to 
identify attractive regions for more detailed study. The figures of merit developed in this paper can also be 
helpful for choosing between available insulation and structural materials as well as providing guidance for 
developing more efficient materials. 
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